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Abstract. A new two-parametric family of mass dist ribution for spherical 
stellar systems is considere d . It generalizes families by iKuzmin fc VeltmannI 
()1972[ ) and by lAn fc Evanj ()2006[) . Steady velocity dispersions are found for 
these models by solving an equation of hydrostatic equilibrium. Axisymmetric 
generalizations of the model are discussed. 



1. Introduction 

In this study we suggest a new family of spherical mass distribution models 
that generalizes models by An Evans (2006, hereafter AE) and mo dels by 
Kuzmin et al. ( Kuzmin &: VeltmannI [T"972l : iKuzmin" fc Malasidzj 119691 1 ■ The 
family depends on two structural pa rameters. It includes Plummer's spheres 
( Plummerl 191 ill . Henon's isochrones ( Henon 19591 ) and the model by Hernquist 
( Hernquistlll990l ) as special cases. 



2. The potential— density pair 

Let us consider the dimensionless potential 
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a — 1 + w{r) 



, w{r) = {l + KrPYlP, K = 0{aP), 
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Here a > 1 and p > are structural parameters. If a = 2, < p < 2 we obtain 
a model by AE, for a > 0, p = 2 we have a model by Kuzmin et al. (1969, 
1972). 

Poisson's equation yields the following expression for density 
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It follows from ([2]) that our models with < p < 2 are cusped (as models by AE 
are) . The density profiles for different values of parameters are shown in Figures 

mm 



The circular speed is found to be 
f ^(r) 



(3) 
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Figure 2. The density profiles of tlie model for various values of a. Note 
that all models are cusped, as p < 1. 
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3. Hydrostatics 

A run of velocity dispersion cr^ can be found from an equation of hydrostatic 
equilibrium 

i^f^ + 2„^?=-^. ,= 1-4. (4) 

g dr r 




Figure 3. Run of velocity dispersion for various values of p. 

The results of calculations for an isotropic velocity distribution (/3 = 0) are 
shown in Figure [3l Central minima will appear in all models with density cusp. 
It can be obtained from ([1]) that 

Then it is possible to find an augmented density g{^) and calculate an isotropic 

distribution function. 

Stability of such models can be studied using the third Antonov law ( Binnev &: Tremaind 



19871 ). namely, if -j-^ > the model is stable against spherical perturbation. 



The validity of this inequality can be established after some laborious calcula- 
tions. 



4. Axisymmetric generalizations 

Using the equipotential method by iKutuzov Sz Ossipkovl ( 198ll ) one can con- 
struct axisymmetric generalizations of the suggested model. We considered a 
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potential of such models 

viR,z) = <P{0, e = f{R,z), (6) 

where ^{S,) is the same function as ([T]) and f{R,z) = const is an equation of 
equipotential sur faces, R, z being cy l indric al coordinates. We considered the 
equipotentials by iMivamoto h Nagail (|l975l ): 



^2 = ^2 _^ ^2 ^ 2 (1 

and bv lSatohl (Il98nl ): 



,)(^Vi2^_£^^ (7) 



^2 = ^ ^2 ^ ^4(l_e)^2+^2 _ (8) 

Here e G [0, 1] is a new structure parameter. For spherical systems e = \. 

We found that for e close to 1 the density is positive for i?, z everywhere. So 
we concluded that such non-spherical model can be used for approximating mass 
distribution in non-spherical star clusters and non-highly flattened galaxies. 




Figure 4. The de nsity profiles of the axisymmetric model with equipoten- 
tials by IMivamoto &: Nagai. ((1975i) for various values of e. 
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